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Abstract. Interactions between dark matter and dark energy, allowing both conformal and 
and disformal couplings, are studied in detail. We discuss the background evolution, anisotropies 
in the cosmic microwave background and large scale structures. One of our main findings is 
that a large conformal coupling is not necessarily disallowed in the presence of a general dis¬ 
formal term. On the other hand, we find that negative disformal couplings very often lead to 
instabilities in the scalar field. Studying the background evolution and linear perturbations 
only, our results show that it is observationally challenging to disentangle disformal from 
purely conformal couplings. 
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1 Introduction 

Observations of the cosmic microwave background radiation (CMB) and large scale structures 
(LSS) have allowed cosmologists to formulate a model of cosmology in which the standard 
model (SM) particles are a subdominant matter form. The model predicts the existence of 
dark matter which only interacts very weakly with itself and the other matter particles. In 
addition, the model requires an energy form with negative pressure, dubbed dark energy, 
which is responsible for the accelerated expansion of the universe in the present epoch (see 
[1] for the 2013 results of the Planck mission). A major task of present day cosmology is to 
illuminate the properties of dark matter and dark energy. 

A simple candidate for dark energy is the cosmological constant. Its biggest drawback 
however, is that its observed magnitude is 10 120 times smaller than the value expected from 
theory, when interpreted as a vacuum energy density. Because of this, cosmologists have 
studied other possibilities, such as dynamical scalar fields, or modified theories of gravity. 
We refer to [2], [3] and [4] for recent reviews. Here, we focus on a union of the two: the 
case of a scalar held as a dark energy candidate, which modifies the force of gravity. In such 
models, couplings to all matter/energy forms are expected unless symmetries exist which 
forbid or suppress interactions, yet, problematically, a scalar held coupled to matter would 
mediate a long range fifth force between the different particles, a force which is not observed 
in nature [5]. Such non-detection implies that the coupling to baryons must be very small, 
whereas constraints on coupling to neutrinos and dark matter are substantially weaker, and 
must be obtained from cosmological observations. And, very recently, evidence has emerged 
to suggest that an interaction between elements of the dark sector is not just plausible, 
but actually favored by current data ([6] and [7]). The analysis was in each case purely 
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phenomenological, assuming a minimal amount of underlying theory, yet it is a progressive 
step toward understanding the nature of these invisible elements of our universe. 

In light of these facts, we dedicate this work to the investigation of dark energy as a 
very light scalar field coupled to dark matter only, and assume all interactions between the 
standard model and the dark sector are negligible. Theories with an interacting dark sector 
have been discussed in the literature extensively, and our ignorance of this sector’s physical 
nature is reflected in the wide variety of interaction types considered; see e.g. [8] - [24] 
and references therein. The setup in which the interactions are only in the dark sector has 
been motivated from theories with extra-dimensions and branes [25]. As such, the couplings 
between dark matter and dark energy are purely geometric in origin in which dark matter 
lives on a brane distinct from the brane on which the standard model particles are confined. 
In many of these works, the gravity sector of the theory is of scalar-tensor form, the scalar 
plays the role of dark energy, and the coupling of the scalar field to dark matter is described 
via a conformal transformation of its geometry - dark matter now responds to curvature 
via an effective Newton’s constant that depends on the local value of the scalar field. The 
interaction is hence termed a conformal coupling. 

As an extension of this idea, we now allow the additional possibility of disformal couplings 
between the two dark elements. Disformal models of gravity, initiated by Bekenstein [26] have 
been attracting much attention recently, particularly with regards to cosmology, see [25] and 
[27] - [35]. These disformal factors have been used in stabilizing scaling solutions in massive 
gravity [36], modifying the speed of gravitational wave propagation during inflation [37], even 
describing electron transport theory in strained graphene [38], and many other ways besides. 
One of the central issues we address in this paper is whether or not cosmological observations 
will allow us to disentangle the effects of conformal and disformal couplings. As we couple 
the scalar field to dark matter only, we avoid the stringent constraints on disformal couplings 
from a host of local tests [39] . 

To demonstrate clearly what we mean by conformal and disformal transformations of 
the dark matter geometry, let us now write down the action for the theory we consider: 

<5 = J { ^ + £(SM) + £ ( de) } + / (1-1) 

where SM corresponds to the visible sector (i.e. the standard model particles) and dark 
energy, CmE)i is described by a quintessence field: 

c {DE ) = -\y a 4>y^-v{ ( j ) ). (i.2) 

The dark matter sector, described by the Lagrangian 

£(DM) = £(DM){9apW), (1-3) 

depends on the metric 

9^ = C{<j))gqb, v . (1.4) 

The functions V, C and D encapsulate our theory’s remaining freedom, that will be specified 
in later sections. C and D go by the names ‘conformal factor’ and ‘disformal factor’ respec¬ 
tively. We see now that dark matter particles follow geodesics determined by g^, and that 
various aspects of these particles, for instance their mass, will now depend on the dark energy 
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field. The functions C and D could also depend on the derivatives of 4>, but we will ignore 
this possibility for simplicity in this paper. 

The theory we have just made concrete is a mathematical realization of generalized 
gravitational interactions within the dark sector. It encompasses a very broad, though not 
exhaustive, number of alternatives to general relativity, which is now just a point in the 
function space: C = 1 and D = 0. The case of D = 0 but C / 1 has been extensively 
discussed, and only recently have cosmologists studied the implications of the disformal term 
as well. In this paper we will be studying models with C / 1 and D ^ 0 and compare them 
to the case of purely conformal ones. In particular, the interaction between the two types is 
investigated numerically, and we will find this interplay has important consequences, namely 
the efficient suppression of one type by the other. 

The paper is organized as follows: In the next section we discuss the evolution of the 
background system, and specify the different choices of free function ( V , C and D) forms and 
parameters used consistently throughout our analysis. In section 3 we will turn our attention 
to the evolution of cosmological perturbations in the presence of disformal and conformal 
couplings, and compute both matter and angular power spectra for various cases. We sum¬ 
marize our findings in section 4. All numerical work, including background simulations and 
both power spectra, is the output of a modified version of the publically available Boltzmann 
code CLASS [40]. Throughout the paper, we will emphasize coupling type discernabiliy; can 
we actually ‘observe’ a purely disformal phenomenon? 


2 Background Cosmology 

This section is split into several parts. Firstly, we write down the background equations. The 
background dynamics are then described in detail. The effective coupling to dark matter and 
the effective equation of state of the dark energy scalar field is discussed subsequently. 


2.1 Equations of motion 

The background spacetime is the standard Friedmann-Robertson-Walker (FRW) metric so¬ 
lution to the Einstein equations for the metric with flat spatial hypersurfaces: 

ds 2 = g l , iU dx} l 'dx lJ = a 2 (r)[— dr 2 + 5ijdx l dx J ']. ( 2 . 1 ) 

Here r is the conformal time and a(r) is the scale factor. For the rest of the paper, dots 
denote derivatives with respect to r. Note that the disformal metric which dark matter 
particles "feel" is given by Eqns. (1.4) and (2.1) as 

ds 2 = g t _ ll/ dx fi dx l/ = C a 2 {t)[—^ 2 dT 2 + 8ijdx l d3?], (2.2) 

where we define a disformal scalar 7 as 


7 : = 



D 

C 




&1V 


(2.3) 


The background value of the scalar field depends only on r. We assume that neutrinos are 
massless in our analysis, and hence the different sectors of the theory are specified by a 
massless relativistic component, r, and baryon component, b, a dark matter component c and 
the scalar field </>. The relativistic species as well as the baryons are assumed to be uncoupled 
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from the scalar and hence the evolution of their energy densities is described by standard 
conservation equations: 


Pr + 3 T~L{px Py) — 0 

Pb + 3PLpb = 0 


(2.4a) 

(2.4b) 


where P r = p r /3. On the other hand, the scalar field obeys the Klein Gordon equation, and 
is now coupled to dark matter via a coupling function Q: 

<j> + 2'H(p + a 2 -yy = a 2 Q , (2.5) 

acp 

where the background form of Q is given by 


Q = 


a 2 C' 


2D(3H<j> + a 2 V’ + + D'(p 2 

2 {a 2 C + D{a 2 p c -j> 2 )) Pc 


( 2 . 6 ) 


The non-conservation of the dark energy-energy momentum tensor implies subsequent non¬ 
conservation of dark matter; energy loss from one species must be mirrored by energy gain 
in the other, and so for the cold dark matter species, c, we obtain 


p c + 3Up c = -Q4> ■ 


(2.7) 


Finally, from Einstein’s equations we present the Friedmann equation, which takes the stan¬ 
dard form: 

o 8nGa 2 , . . _. 

H = - 2 - (pr + Pb + Pc + Pde) j (2-8) 

where % = a/a and p& e = (j) 2 /2a 2 + V{4>). 


2.2 Analysis of the dynamics 

The dark energy - dark matter interaction encoded in the equations above describes a peculiar 
scenario. For the dark energy field, the coupling contributes to an effective potential which 
depends in general on (j) an d </>• The dark matter then gains and loses energy as the geometry 
described by g^ is stretched and distorted; the conformal factor dilutes the dark matter over 
space-time by modifying the isotropic expansion it feels, C 1 ^ 2 a, while the disformal factor 
distorts dark matter particles’ light cones. The nature of this energy transfer process will 
depend crucially on how we specify our three free functions C((f>) and D (</>). Suffice 

to say for now that we require our cosmology to be empirically plausible - observation tells 
us our dark energy field must roughly resemble a cosmological constant, and the scalar field 
must evolve very slowly to account for the accelerated expansion. To be specific, in this paper 
we study the following forms for V. C and D: 


■o- 

SCL 

03 

II 

(2.9a) 

C = C Q eP c *, 

(2.9b) 

D = M D A e pn4> . 

(2.9c) 


This choice represents the simplest extension to the standard coupled quintessence scenario, 
for which C and V are given as above and D = 0. 
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# 

Name 

Pc 

Pd 

M d 

x behavior 

1 

Uncoupled 

0 

0 

oo 

x = 1 V r 

2 

Conformal 

-0.2 

0 

oo 

decreasing with r 

3 

Disformal 

0 

0 

My 

increasing with r 

4 

Mixed 

-2 

2 

My 

single stationary minimum 


Table 1. Description of the four models used throughout this paper. The model parameters are 
defined in Eqns. (2.9). The function x is defined in Eqn. (2.11), and ‘Name’ corresponds to the 
model’s label in plot legends. Note that in the fifth column, Mp = oo simply represents the limit for 
which the disformal coupling vanishes. For a more detailed discussion of the last column, we refer to 
Sec. 2.3. 


We can, without loss of generality, set Co = 1, as this parameter simply corresponds to 
a global redefinition of units. Such a choice does not affect the dynamics. The dark energy 
scale, My, is taken to be a fitting parameter that must be tuned such that our final time 
boundary conditions agree with measurement of the universe today, that is to say fl^ = 0.68, 
fl c = 0.27 and fib = 0.049. Typically we find My ~ Ho ~ nreV. The conformal coupling 
is dimensionless, but the disformal factor introduces a new scale into the model. The case 
M^ 1 —> 0 corresponds to the standard coupled quintessence scenario, and the opposite limit 
where Md —> 0, it turns out, is actually an uncoupled limit, regardless of the form of C. This 
unexpected feature is a consequence of a suppression effect to be clarified shortly. In between 
these limits, we find disformal effects leave an observable imprint on cosmological observables 
that is maximal if Mo ~ My. We study a variety of models with different values for Md, (3q 
and Pd, and the parameter combinations we consider are summarized in table 1, where the 
meaning of the final column we will specify in the next section. 

In Fig. 1 we show the evolution of the coupling function, in which the most prominent 
feature is an early time suppression induced by the disformal factor. The coupling effectively 
‘switches on’ during some past epoch, quite late in the universe’s evolution, and before this 
time it is in fact completely negligible. We see the effect this has on the evolution of the dark 
sector’s energy densities in the next figure. Compared to the purely conformal case, the scalar 
field receives no great kick at early times when a disformal factor is included, and though 
/3c is here a factor of 2 larger than current experimental upper bounds, the field mimics a 
cosmological constant throughout the majority of this universe’s simulated lifetime. 

Probing different free functions and parameters we find this suppression to be no lucky 
coincidence of the mixed model, but seems a general property of a disformal factor included 
in almost any cosmology. We can see why this must be the case by examining Eqn. (2.6). 
A non-zero D means the presence of a term proportional to p c appears in the denominator. 
As long as the disformal scale Md is of the same order as the dark energy scale or less, this 
term will continue to make the coupling negligible until dark matter is roughly of that scale, 
i.e. today. We will show in later sections that the same screening effect holds too for linear 
perturbations, for the very same reason. 

We have arrived at the first main result of the paper: in a cosmological setting, a 
disformal factor can suppress a conformal contribution at early times. The key point is 
that, for certain orders of magnitude of the disformal factor, the value of Q and its linear 
perturbation, 5Q (whose exact form is given in section 3), are very much diminished for 
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Figure 1. The scalar coupling function, Q , for background solutions of various models, plot against 
the scale factor a where a o is the value of a today. The free functions are defined in Eqns. (2.9), and 
/ 3y = —2 for all curves. The rest of the model parameters are given in table 1. 



Figure 2. Evolution of energy densities of: dark matter (dashed lines) and dark energy (solid lines) 
for background solutions of various models, plot against the scale factor a where ao is the value of 
a today. The free functions are defined in Eqns. (2.9), and /3y = —2 for all curves. The rest of the 
model parameters are given in table 1. Here the dark matter curves for both models coincide. 
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the majority of the universe’s evolution - in the presence of disformal couplings, significant 
conformal ones are not necessarily in disagreement with cosmological observations. 

2.3 The characteristic coupling function, x 

The expression of Q is not simple to analyze, and collecting all our coupling effects under the 
obfuscated umbrella Q has somewhat obscured the physics. It’s form is necessary for finding 
numerical solutions, but for the analytics we can do better. To elucidate the effects conformal 
and disformal couplings have on cold dark matter, we will now define a new variable x that 
greatly simplifies the analysis, and, as it turns out, the dark matter equation (2.7) will become 
easily solvable. In fact, this remains true for any species whose exact solution can be found 
in ACDM, for example photons. We relegate the details and general case to appendix A, but 
for pressure-less dark matter we obtain 


Pc = Po, c a 3 x , 

where we have defined the quantity 

C^/l 

x := —T72- ‘ 

Co ho 

Another useful quantity will be the derivative of x: 

x Q4> 

x p c 

which we express in terms of a rate. Then, looking again at Eqn. (2.7), we see now the 
evolution of p c as a competition of rates: that of the Hubble expansion rate, Ti, and the rate 
of interaction with the scalar field, x/x. 

The positivity of x follows naturally from the fact that both C and 7 separately must be 
positive. This condition is defined by the metric (2.2); we must preserve causality, or suffer 
the consequences. Throughout the course of this paper, it will become clear that the whole 
system can be characterized by x and its derivative - at least, in terms of observables - at 
both the background and perturbative level. We can already see this to be true at the zeroth 
order. As x contains only background quantities however, it is quite remarkable that this 
remains true at first order. 

With x defined, it is now time for us to address the issue of our large free function space. 
What we are looking for in this study is general characteristics of conformal and disformally 
coupled dark matter, not idiosyncrasies corresponding to specific choices of the functional 
forms and parameters of V. C and D. To make this step toward more comprehensive conclu¬ 
sions, we first notice that - as we have already stated - both C and D do not actually work 
independently, but affect the system jointly through x. We will then partition observationally 
distinct models based on the behavior of their respective x function in conformal time. 

As previously stated, we keep the models presented here realistic, with observables like 
the CMB anisotropies close to their measured values, and so we must work in the slow roll 
regime. This means that not only should V be a relatively shallow function, but so too C, as 
it is also able to drive the field. The disformal factor, however, induces a damping in the field 
dynamics, and we find it can not push the field by itself, but rather hinders its movement. 
This damping will give us some more leeway in how shallow C and V can be. What this will 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 
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ultimately mean is that we do not consider scaling solutions or attractors; this work is not 
aimed at solving the coincidence problem, rather, we find an alternative notion of naturalness 
is manifest here: the general inclusion of a disformal factor serves to push an arbitrary coupled 
cosmology toward one with a cosmological constant. 

Given what has just been said, we suggest that for a qualitative first study, it will be 
enough to consider just four distinct models: 

1 . uncoupled quintessence, where x = 1 V r, 

2 . x is a decreasing function of time, 

3. x an increasing function of time, 

4. x has a single minimum. 

We can now comment on the final column of table 1 mentioned earlier. In particular, it is 
not of vital importance that V is of exponential form, a power law e.t.c. What does matter 
is whether x is pushed upward, to larger values, or made to roll down to lower ones. So, only 
the direction of the slope of C relative to V affects the evolution of cold dark matter. 

We show the x behavior for the four models in figure 3. Note again the defining feature 
of the disformal term is that x - and hence the coupling - vanishes for early times, while 
for the conformal model it rapidly diverges. At early times then the distinction between 
conformal and disformal effects is strikingly clear, but at late times however, this is not so. 
We have used here a conformal only model to produce a decreasing x function model, but we 
could have achieved this by other means. For example, were we to pick a more complicated 
choice of D function, and a different potential, we would get qualitatively the same late time 
behavior. At least at the background level, it is this late time x behavior that is observable, 
as long as dark energy remains sub-dominant, and we will see in the next section why we 
have categorized our models based on this criteria. 

2.4 An effective equation of state 

The dynamics of our gravitationally coupled system are in general quite complex. There is 
energy transfer between the elements of the dark sector that depends not only on the dark 
energy held, but also it’s first and second derivatives. We anticipate that when interpreting 
data, cosmologists will use a much simpler parameterization. This assumed model is most 
often of a non-interacting dark sector, where dark matter is pressureless dust and dark energy 
some fluid with an open equation of state. Following [41] we now reformulate our theory at the 
level of the zeroth order equations of motion to fit this neat picture, and define an effective, 
or apparent, dark energy equation of state uyj. 

In the Friedman equation we first perform an effective splitting between the two dark 


components: 

„ , 2 87 t Ga 2 , o s \ 

n = ( p b ,oa ° T p c ,oa ° T Pde.eff ) 

(2.13) 

with 


—3 

Pde,eff P(j) Pc,0^ “1“ Pc* 

(2.14) 


Then, taking the time derivative of /9de,eff and defining 


Ale.eff 474 (IT w: e ff) , 
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a/a 0 


Figure 3. Evolution of the characteristic x function for the four models described in table 1, with 
Pv = —2 for all curves. The model parameters are defined in Eqns. (2.9). 


we find 

iu e ff = ^ with ptj, = X — V, (2-16) 

Pde,eff 

where X = (f/ 2 /2a? is the kinetic energy of the field. This quantity is the apparent equation 
of state of dark energy an observer would infer, if the assumption is made that the energy 
density of dark matter scales scales like a~ 3 . In our theory of course this is not the case, and 
ic e ff / w$ in general. 

What can we expect to observe in the behavior of this new effective system? An inter¬ 
esting first question to ask is: will we see phantom behavior? Using Eq. (2.10), it is simple 
for one to derive the following phantom condition: 

W e ff < —1 2X < /9 C) oO _3 (l — x ) (2-17) 

where we recall that x is normalized to unity today. Its clear that the evolution of x will 
dictate whether or not we see the effective dark energy cross the phantom line, and this is 
directly related to the coupling of the underlying true model: if x is an decreasing function 
phantom behavior is impossible, and in this scenario energy flow is from dark energy to 
dark matter. Conversely, energy flow in the opposing direction (x is an increasing function) 
will propel the universe toward even greater expansion, as the relative contribution to the 
cosmic inventory from the vacuum energy will grow. Clearly, this system should not exhibit 
the standard instabilities expected from true phantom dark energy models, as it is simply a 
phenomenological re-parameterization. 

For our four chosen models in table 1, we display the effective equation of state as 
discussed above in the redshift range accessible to the proposed Euclid satellite and the 
majority of current redshift galaxy surveys. While the disformal model here goes divergently 
phantom, for the others, the phantom line is never crossed. In the pure conformal case, 
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Figure 4. Evolution of the apparent equation of state, w e s , defined in Eqn (2.16) for the four models 
described in table 1. The model parameters are defined in Eqns. (2.9). 


the effective system tends further from ACDM, toward the boundary between acceleration 
and staticity. Das et al. [41] however find a conformal model that replicates our disformal 
model’s divergently phantom behavior, and so the take home message is emphatically not that 
disformal couplings induce effective phantom behavior while conformal ones do not. Rather 
the point is that if x is increasing with time, phantom behavior will likely ensue - a direct 
result of the phantom condition (2.17). 

At the background level then, the coupling has a nice interpretation as a variable dark 
energy equation of state. This correspondence is best illustrated through its effect on lumi¬ 
nosity distances, cIl■ In figure 5 we show the luminosity distance difference ratio for our four 
models, defined as: 


A d L 
dh 



dh,i 


dh 


,uncoupled 


dh 


, uncoupled 


(2.18) 


As our intuition suggests, energy flow into the scalar field accelerates expansion, causing 
observed objects such as supernova at a fixed redshift to appear further from us than for the 
uncoupled case. 

To conclude this section, we have shown that conformal and disformal effects can not 
always be distinguished when only dark matter is coupled. Whether they can or not typically 
depends on the epoch in question: at early times the distinction is clear, as disformal contri¬ 
butions in general suppress conformal ones; at late times the two act together through the 
x function, and whether energy flow is into or out of dark matter tells us nothing about the 
underlying behavior of C and D, nor will any observed phantom behavior. What defines early 
and late times in this context is the new scale introduced by the disformal factor, Mo. When 
the dark matter energy density becomes comparable to that scale, the coupling switches on 
and begins to influence the field. 
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Figure 5. Evolution with redshift of the fractional difference in luminosity distance, Eqn. (2.18), 
between the four models described in table 1 and the uncoupled case. The model parameters are 
defined in Eqns. (2.9). 


The reason why the two types of coupling become indistinguishable at late times is 
ultimately because dark matter is pressureless. What sets disformal factors apart, we have 
seen, is that they warp light cones and shift causal structure, but this has little effect when 
the particle constituents of the coupled fluid is cold dark matter. 

3 Evolution of cosmological perturbations 

We now turn our attention to the evolution of cosmological perturbations in our theory. We 
begin by first writing down the perturbation equations and subsequently discuss predictions 
for cosmological observables, such as the CMB anisotropies and matter power spectrum. 
Along the way we will try to be categorical about the various effects induced by the couplings 
for the different models; will there be anything about these spectra that is characteristically 
disformal? 

To be concrete, we will be working in the Newtonian gauge. To avoid confusion, we 
reserve 6 to denote matter density contrast: S := 5P the pressure perturbation and 6(f) 
is the perturbation of the scalar held. We denote by 5 a general perturbation operator. The 
perturbed Einstein frame line element in the chosen gauge is: 



(3.1) 


which means that dark matter particles follow geodesics described by the following perturbed 
space-time: 


ds 2 = Ca 2 (r)[—(1 + 2A)'y 2 dr 2 + 2(diB)^dTdx l + (1 — 2 E)5ijdx l 'dx^] (3-2) 
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where A. B, and E are functions of the dark energy field background and perturbation values. 
Their exact forms are: 


—I 


(3.3a) 

(3.3b) 


(3.3c) 


Its now clear that only disformal factors induce off diagonal perturbations in the metric which, 
we will see, affect the velocity field perturbations. 

For the visible sector we again neglect that neutrinos have mass. The perturbation equa¬ 
tions of relativistic, r, and baryonic, b, matter follow from the energy-momentum conservation 
equations and are given by 


5 r 

9 r 


fi b 

d b 


(1 + W r )(0 r — 3$) — 377 ^ ~ W r^j fir 

(3.4a) 

77(1 3w r )0 r UV 6 r + kH + 5Pr ^ Pr k 2 5 r 

1 + VJ r 1 + W r 

(3.4b) 

0 b + 3T 

(3.4c) 

77 e b + fc 2 <h. 

(3.4d) 


Perturbations in the dark energy field, 5(f), evolve according to the perturbed Klein Gordon 


equation 

5j> + 2 U54> + ( k 2 + a 2 V")5(f> = d>(T + 3$) - 2 a 2 (V' - Q)T + a 2 5Q (3.5) 

and perturbation of Q is given by the cumbersome expression [31] 

5Q = - nf C 2 -+ B26 + B ^ + B 4^ + (3-6) 

+ UyCL^p c — (p J 

where 

51 = ^ - ZDHi - Da 2 (V' - Q) - Dft , (3.7a) 

5 2 = 3£ty, (3.7b) 

B3=6DW0 + 2^(g-g + k), (3.7c) 

Bi = -3DH - 2D4, (§ - g + £) , (3.7d) 


+ (a 2 C' + Tl'a 2 ^ - T»> 2 ) — . (3.7e) 

Pc 
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The non-conservation of dark matter induces, at the perturbative level, factors of Q and 5Q, 
in its conservation equations, which become 


4 

4 


Q 

Pc 


cj)5 c — —SQ — —6(j), 

(3.8a) 

Pc Pc 

- —j>0 c - —k 2 54>. 

Pc Pc 

(3.8b) 


An important point has to be clarified before we continue: are we justified in setting the 
Einstein frame pressure and shear of the coupled dark matter to zero? After all, its been shown 
previously that disformal couplings can transform the equation of state of a coupled species 
to one dependent on the field, <f> (see [42], [43] and appendix A). In the appendix we present 
the full set of frame transformations between Jordan and Einstein frame matter variables 
(or equivalently, between uncoupled and coupled variables respectively) 1 . In particular, for a 
general coupled species: 


5P 


1 

C^f 


5P 



2 SC <5 7 \ 

- + j) 



(3.9a) 

(3.9b) 


So, as dark matter has vanishing pressure in the frame for which it is uncoupled, we may 
indeed set SP C , II[? = 0. 

The perturbed Einstein equations in our theory take the form as in standard ACDM. 
We don’t quote them all here, but instead present just those important for our analysis 
of structure growth. In particular, we will make use of the 00—component of the Einstein 
equations 

k 2 $ + 3 n^ + = -4vr Ga 2 5p (3.10) 

and the i ^ j— component, which leads to 


4>-T = 0 V C, D, (3.11) 

as we ignore anisotropic stress (we remind the reader that we ignore neutrino masses in our 
analysis). Our theory thus predicts no gravitational slip (?/ := <h/T = 1), independent of the 
coupling. 

With the perturbed equations written down, we first ask whether the disformal factor 
suppresses a conformal one at the level of linear perturbations. Before getting to the figures, 
we can already guess that it will be the case from Eqn. (3.6); looking at the denominator we 
see the same Da 2 p c term that was responsible for the suppression at zeroth order. Still, its 
important to be concrete, and so we demonstrate this suppression in Fig. 6 and Fig. 7. 

When we compare figures 1 and 6 a similarity immediately leaps out: one plot appears 
to resemble the negative of the other. Additionally, the curves for 5Q appear to rise and fall 
gradually, over background timescales, rather than, for example, the fast oscillations in the 
field perturbations (figure 7). This simplicity may come as a surprise when juxtaposed with 
the arduous complexity of the 5Q equation from which the curves have sprung. In fact there 

1 By Jordan frame, we mean the frame in which the action (1.1) is rewritten fully in terms of g M „, whereas 

the Einstein frame is specified by (1.1). 
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Figure 6. Dependence on the scale factor of the linear perturbation in the coupling function, 5Q, 
given by Eqn. (3.6). The three models shown are described in table 1, with equations (2.9). 



Figure 7. Evolution with the scale factor of density contrast absolute values in: the dark matter 
energy density (dashed lines), and dark energy density (solid lines). These curves correspond to modes 
with wavenumber k = 0.3 Mpc -1 . The models shown are described in table 1, with Eqns. (2.9). 


is a strong relationship between the background and perturbed coupling function Q , and the 
damped oscillatory behavior of 5Q at a/ao ~ 2 X 10 -2 (figure 8) betrays a important aspect 
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Figure 8. We show here, for models only including a disformal term, how the system approaches 
Eqn. (3.12). Dashed lines: QS C , where Q is defined by Eqns. (2.6), and S c is the dark matter density 
contrast. Solid lines: SQ, given by (3.6). The models shown are described in table 1, along with Eqns. 
(2.9). 


of the coupling - the system at late times is drawn to a solution where [30] 

SQ ~ QS C , (3.12) 

which we identify with the limit in which all dark energy perturbations in SQ vanish. The 
approximate equivalence, Eqn. (3.12), does not hold at very early times but becomes progres¬ 
sively more accurate later on. In Fig. 9 we plot the oscillating part of SQ for an exclusively 
disformal model. We see clearly that both the oscillation and damping time scales depend on 
the mode’s wavenumber, k, which we expect, and the new scale introduced by the disformal 
factor, Mjj. In the conformal limit, characterized by Mjj —> oo, we find in general they 
disappear as their period tends also to oo: a conformal factor in general introduces no scale, 
so a purely conformally coupled theory does not exhibit these oscillations. 

To briefly summarize by way of analogy, at the background level the new disformal scale 
determined an epoch in which the coupling is effectively ‘turned on’; at the perturbation level, 
Mfj now sets a damping and oscillating timescale for perturbations in the dark energy held. 
The introduction of this scale, we can now see, is primarily what sets apart the effects of 
conformal and disformal factors on the scalar held dynamics at the zeroth and linear order 
levels of the theory. 

3.1 Growth of large scale structure 

We have seen some evidence to support that, at least in principle, conformal and disformal 
effects can be separated. The new important feature is of course the new scale, Mo. We 
now turn to the pressing question of observables: does the new scale also leave recognizable 
imprints on the formation and subsequent growth of structure? 
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Figure 9. Oscillations induced by a disformal term in the perturbed coupling function (3.6) for models 
where /3y = —2, Pc = 0, Pd = 0 and Md = Mv/ T|!/ 4 . Terms in 5Q proportional to dark matter 
perturbations only, Q5 C , have been subtracted to isolate the dark energy perturbation fluctuations. 


The key quantity we are interested in is the growth factor, f(z), which is a convenient 
parameterization of linear growth. In the literature, the growth rate is defined as: 


f(z) 


d In _ <5 m 
d In a Udm ’ 


(3.13) 


where we notice the definition is with respect to all matter, 5 m := 5p m /p m for p m = p c + p b , 
not just dark matter. As our theory only couples dark matter to the scalar, we are now faced 
with an interesting question: how will the composite fluid of dark matter and baryons cluster 
into structure, if both species feel different (effective) gravitational forces in general? We 
would like to calculate the growth equation for not just cold dark matter, but for all matter. 
It is this total matter quantity that influences the gravitational potentials, which lens distant 
galaxy and CMB light. As a first step then, we now define a ‘baryon bias’ parameter b, as: 


5 b = b6 m (3.14) 

where b will in general depend on time (and possibly scale, but since the scalar field is nearly 
massless, we find that b in our model does not depend on the wave number k ; the situation 
would be different for different potentials which are not of quintessence form). 

If we plot the late time evolution of this bias (figure 10) we see as we might expect that 
dark matter and baryons cluster at different rates, reflecting the underlying variations in each 
species’ experience of gravity. What, primarily, we glean from the plots is that to reasonable 
accuracy, 5 b — S m ; from this we directly infer also that 5 C — 6 m . We will return to address 
the validity of this assumption very soon, but for now, using this, along with sub-horizon 
approximations and Eqn. (3.12) let us derive the linearized growth equation for matter when 
cold dark matter is gravitationally coupled to dark energy. 
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Figure 10. How the late time bias, as defined in Eqn. (3.14), varies with redshift, z, for the models 
described in table 1. The free functions that define the parameters displayed in the table are given 
by Eqns. (2.9). 


Let us take for granted momentarily that 5 C — S m . By the additivity of the stress energy 
tensors and baryon conservation, the evolution of all-matter perturbations (dark matter + 
baryons) can be derived from: 

s (V^TD = <5 (yx u + V^r) = 8 (QP) (3.15) 

which gives the growth equation: 

8m T T~L e f[8 m = 47rG e ffa Pm8m (3.16a) 


where we have defined: 


H cS :=H + i^ 

X Pm 

G eff ■.= G+-X- (-^\ 

47T(/> 2 \XPmJ 


= G + 


L9l 

4tt p 2 m 


(3.16b) 

(3.16c) 


which is valid for k 7~L. We note that while the error I Sb I is of the order ~ 1 — 3% for the 

Om 

models considered, the error this induces in the growth equation (3.16a) turns out to be only 
ever as large as ~ 0.1% in general, and usually substantially smaller. The error propagates 
through the derivation in a favorable way, affording us valuable comparison between the true 
evolution of h m and our simplified growth equation to an accuracy sufficient for this study. 
We stress that, of course, evolution of cosmological perturbations is inextricably linked to 
evolution of the background. 

We are now in a position to examine the growth rate, (3.13) for our various models. 
First to note is that both T~L e g and G e g contain only background quantities that have no k 
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Figure 11. The fractional difference in the growth rate f(z), Eqn. (3.13), for the models of table 
1, against redshift, z. The difference is with respect to the uncoupled case, i.e. A/// := (fi — 
/uncoupled)//uncoupled where i runs over the four models. 


dependence - any departure from general relativity here will be scale invariant. This is due 
to our choice of potential for the scalar field, which is of quintessence form and the field 
is nearly massless. With this choice, our theory predicts that, like ACDM, measuring the 
growth at different length scales (within the quasi-static regime of course) will not lay bare 
the novel features of the coupling presented here. If the oscillations depicted in Fig. 9 were 
perhaps to have survived till today, this would change the story: a major observational test to 
distinguish disformal couplings would then be to see how these fluctuations depend on scale, 
k , and thus probe the value of Mjj itself. In general though this is not the case; the severe 
damping present in all models we have considered show that sustaining these oscillations long 
enough to observe them is difficult to achieve in practice, and so highly unlikely in reality. 
What looked in previous sections like a tool to measure disformal couplings turned out to be 
just an fleeting fluctuation, completely intractable empirically. 

Looking at the curve for the purely conformal case in Fig. 11, it seems that, because 
both % and % e g are suppressed by the coupling (the background expansion rate is slowed as 
energy is transferred from dark energy to dark matter, exemplified by the effective equation 
of state for this model) growth is enhanced by the coupling. For the disformal only model, 
cosmic expansion and the extra friction felt by 5 m (% e ff) is enhanced but we see that the 
growth rate is largest for the purely disformal model. 

In General Relativity, the growth rate defined in Eqn. (3.13) is simply related to the 
growth of the gravitational potential <h. This relation is slightly modified in the coupled 
quintessence scenario, as we will now show. We begin by using Eqn. (3.10) in the quasi-static 
regime (valid in the sub-horizon limit, deep inside the matter dominated epoch): 

A: 2 <3? = —AirGa 2 ( 5pb + 5p c + 5p^ e ). (3-17) 
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Figure 12. Evolution of the growth rate /$, as defined in Eqn. (3.19), as a function of redshift, z, 
for the models of table 1. 


It turns out that the dark energy perturbation is negligible compared to the contributions 
from the baryons and dark matter (we have checked this numerally). Then, remembering that 
5p m = Spb + 5p c and using the background equations, we can derive the following equation 
for $: 

$ _ n, , <*rn . X p c 

— — + -- 1 -. 

$ Om X p m 

Here, x is the quantity defined in Eqn. (2.11). If we define 

f _ dln(a$) _ (a<3?)' 

U := din a = 

then we find 

p c din x 
/$ — / H wv . 

p m d m a 

We see that in the uncoupled case, for which x = const = 1, /$ and / coincide. We plot 
the behavior of /$ in Fig. (12). Whereas / characterizes the growth of the density contrast 
in matter, /$ characterizes the growth in the gravitational potential <h. For /$ now, the 
departure from the uncoupled case is significant, especially at very late times - well within 
the reach of current redshift galaxy surveys. Crucially, this new growth rate gives us a direct 
measure of the gravitational potential $ which describes the shapes of gravity wells into which 
galaxies fall, and the CMB is lensed. 

Having diverged slightly, we now return to our central question: are conformal and 
disformal effects at all separable observationally at the perturbation level? Unfortunately, it 
would appear not. Just as for the effective equation of state, it seems that what x is doing at 
very late times (when the coupling suppression has already been lifted) is what dictates how 
either growth rate, / or /$, will behave: the disformal oscillations, characteristic of a newly 


(3.18) 

(3.19) 

(3.20) 
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Figure 13. Top: angular power spectrum, Cp against angular scale, /, for the four models described 
in table 1 with parameters defined in equations (2.9). Bottom: fractional difference between the 
various models and the uncoupled case, i.e. ACj/C* := (Cj.i - Cz iUnc0 upied)/Cz,uncoupled- 


introduced scale, die out long before the present day, along with any hopes of discerning 
between the conformal and disformal factors. Again, just as at the background level, the 
dynamics of late time growth are determined by the late time behavior of the function x - 
a degeneracy between C and D - no matter how it comes about. Whether the two coupling 
types can be separated, again, depends on the epoch. Earlier on, a distinction is manifest, 
but at the later stages of universe evolution, the two act together, and the distinction blurs. 
As before, ‘later’ is defined by the scale Mp relative to the evolving Hubble scale. 

3.2 The power spectra 

We will now discuss the predicted CMB anisotropies and the matter power spectra. In 
Fig. 13 we show the angular power spectrum for the CMB anisotropies. For the type of 
models discussed in this paper, if a particular model has a reduction in power at low l values 
it will have an enhancement at large ones, and vice versa, with respect to the uncoupled 
case. At small multipole values the angular power spectrum is dominated by the integrated 
SachsWVolfe (ISW) effect, which is dependent on the late time behavior of the large scale 
gravitational potential. If a particular model undergoes enhanced expansion at latest times, 
quantified in an earlier section by an effective dark energy equation of state crossing the 
phantom line, then the gravitational potential on large scales decays. The corresponding 
low l anisotropies in the CMB are reduced as a consequence. For the same model the opposite 
happens for large multipoles: the anisotropies are enhanced. Since we have assumed fixed 
the boundary conditions for the cosmological parameter (fide,o = 0.7, 77o, etc) at the present 
time, the cosmological parameters differ (slightly) for the individual models at the time of 
decoupling. This results in different relative heights for the peaks at high multipoles. 
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Figure 14. Top: matter power spectrum, P, against mode scale, k, for the four models described in 
table 1 with parameters defined in equations (2.9). Bottom: fractional difference between the various 
models and the uncoupled case, i.e. A P/P := (Pj — P, mC oupied)/-P-uncoupied • h is the dimensionless 
reduced Hubble constant, defined as Hq = 100 hMpc -1 .. 


The matter power spectra are shown in Fig. 14. For all models discussed here, there is 
an enhancement of power on small scales (large wave numbers), which is a direct result of the 
new scalar interaction between dark matter particles, which is always attractive. The peaks 
of the baryonic acoustic oscillations are shifted, which is due to fact that the cosmological 
parameters are different at time of decoupling. The couplings do not directly influence the 
position of the peaks. Strikingly, as for the Cf's. the mixed model lies very close to the 
uncoupled case. As we found before, the conformal contribution to the effective coupling is 
suppressed by the disformal factor. 


3.3 Disformal Instabilities 

Up to this point we have studied cases for which D > 0. We will now discuss models in which 
D is negative. As we will see, we find runaway growth of perturbations in the scalar held, 
which we feel justify our choice to neglect this case for the entire preceding study. To see how 
a disformal theory can become unstable in its perturbations, we can first re-write Eqn. (3.5) 
in a more suggestive form, for simplicity treating the case where C = 1, and D is constant: 


5(f) + 2T-L 


i-ii+^u 

2 a/a 


6<j) + (r + a z V")[ 1 - £]5<f> = S(a, k ) 


where, using that |1 + ( p c D ) x | 3> | |, we have defined the parameter £ as 


£ ; = 


1 


1 + 


l ’ 

PcD 


(3.21) 


(3.22) 
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Figure 15. The instability function defined in (3.22) for two purely disformal models. For both 
curves, fiy = —2, /3c = 0, Pd = 0 and Mjj = Myj rj/ 4 . 


and all terms that do not contain 5(f) or its derivatives are collected in the source term S(a, k). 
The homogeneous solution evolves according to 


H a + 




Ha + UJ 2 5(t> a = 0 , 


(3.23) 


with uJ 2 := ( k 2 + a 2 V")[ 1 — £]. We see that uj 2 is always positive if D is positive, and will 
become negative if D is negative. Eqn. (3.23) and Fig. 15 demonstrate this clearly: for 
positive D, £ < 1 and the effective oscillator frequency uj is real; for negative D, £ > 1 and 
the frequency can become imaginary - the <5</> a solution becomes an exponentially growing 
function. So much for perturbations, but the background system is also unstable here. An 
epoch where ( p c D)~ l ~ —1, and hence the system traverses a pole, we can see is almost 
guaranteed to occur for negative D at some time r, whether this happens at higher redshift 
or in our future. 

The above analysis pertained to a simplified model, however, what we have just shown 
encourages us to conclude that such an instability can occur wherever the disformal factor 
goes negative, and this is indeed what we find numerically. Our results suggest that theories 
with negative disformal factors can not be considered as viable. 


4 Conclusions 

In this paper we have studied the observational consequences of an extension to the coupled 
quintessence scenario, incorporating disformal terms. By keeping dark energy in the slow roll 
regime, we have focused not on solving the cosmological coincidence problem, but rather on 
searching for observable signatures of realistic (near ACDM) coupled dark sector theories. 
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Studies like these are an imperative when so little about the dark sector is known and so little 
can be assumed. 

An important result of our study is that tension between a model with large conformal 
coupling (/3c in our notation) and data can be alleviated by the addition of a disformal 
interaction. This is because the disformal contribution very effectively suppresses the coupling 
function Q and it’s linear perturbation SQ for a significant portion of the universe’s lifetime. 
The suppression is clearly manifest in both power spectra, as the predictions of the mixed 
model are very close to those of ACDM, although /3c is of order one. We also find that, 
reformulating the theory in terms of a disformal scalar, the conservation equation for a coupled 
matter species is solvable, and we have used the solution to derive a condition indicating 
whether or not a given coupled dark energy model could be interpreted as exhibiting phantom 
behaviour. 

Our analysis of the perturbations tells us that, as in the standard coupled quintessence 
model, the disformal term does not affect the gravitational slip. Additionally the growth 
rate of the matter density contrast and the growth rate of the gravitational potential no 
longer coincide in models with couplings - we expect this will provide an observational key 
to breaking degeneracies between information contained in the CMB, gravitational lensing 
and LSS. We furthermore find that a negative disformal coupling D generically induces dark 
energy instabilities: perturbations in the scalar held will eventually grow quasi-exponentially. 

The results of this paper suggest that it is very difficult to discern between conformal 
and disformal effects when using only background observables and first order cosmological 
perturbations (e.g. CMB anisotropies or the matter power spectrum); an analysis of the 
non-linear regime on small scales it seems will be necessary to look for a ‘disformal smoking 
gun’. 

This preliminary study, far from complete, still must make direct contact with data. We 
leave this task for future work, where we intend to use CLASS’S Monte Python to confront 
the expansive data set open to cosmologists today, ft would also be interesting to investigate 
in how far the mixed model discussed in this paper can be mimicked by a time-varying 
conformal coupling, such as those studied in [44], Further, more of the theory’s functional 
freedom (£(de)> C, and D) must be explored. How robust will our conclusions be under 
relocations within this function space? In addition, the question of the impact of quantum 
corrections needs to be addressed. Models with both conformal and disformal couplings 
between dark matter and dark energy are motivated from string theory (see [25]), and in 
these models, the scalar held is a DB1 held where the functions C and D are specihed by the 
extra-dimensional space. As a consequence, the effective coupling Q has a different form and 
behaviour. We will turn our attention to such models in future work. 


A Frame transformations: background 


We here derive the set of transformations relating background matter variables in the Jordan 
and Einstein frame of a generally coupled theory (see also [42] and for a discussion in con¬ 
siderable detail including perturbations see [43]). This will allow us to solve the background 
conservation equation for a coupled species in the Einstein frame. We begin with the defini¬ 
tions of the stress energy momentum tensors, and from there, using the map between them, 
compute the transformation rules. In the Jordan frame: 


2 KV—g£(DM)) 
V^g $9nv 


(A.la) 
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for which we can now impose a perfect fluid description, hence defining a Jordan frame energy 
density, p, velocity held u ,l : and pressure, P: 

T' w = (p + P)u»u u + Pg,j V . (A.lb) 

In the Einstein frame: 

T nv , = 2 5 (V~9^(DM)) 

V~9 S 9»v 

where we then define: 

= (p + P)u»u u + Pg 

A map between the two objects can readily be derived [30]: 

dfja/3 jia/3 
^9^v 

= C^l + ^^T^ 

= C 3 7 r^, (A. 3) 

where we recognize the disformal scalar 7 from Eqn. (2.3) which parameterizes the relative 
contribution of the disformal factor. Note as D —> 0, 7 —> 1. Now, choosing the Einstein 
frame line element to be of Friedman Robertson Walker form: 



(A.2a) 
(A.2b) 


ds 2 = g^dx^dx 1 ' = a 2 (r)[— dr 2 + 5ijdx l dx J ], (A.4) 

means that, using Eqn. (1.4), we get, in terms of the disformal scalar 

ds 2 = g^udx^’dx 1 ' = Ca 2 (r)[—q 2 dT 2 + 8ijdx l dhP\. (A.5) 


Given the background metric choice, the fluids are homogeneous, hence no forces are 
exerted between elements of the fluid - each element follows a geodesic dictated by the metric. 
This means the 4-velocity held can be computed directly from it: = dx^/\ds\ and u^ = 

dx^/\ds\. Using this and the map (A. 3) we get the full list of variable transformations between 
the Jordan and Einstein frame background quantities: 


where 0 = 3-^- 

a z 


© 

p 

P 

w 


1 


C 1 / 2 7 

1 

C^i 

7 

C 2P 

_^p 

C 2 7 

1 


u ** = 


1 


tW 


7 


0 + 


C 1 / 2r ya 

3(7,0 V’ 
2 C 


% 


3 H, and p := p2X = <$>ja. 


(A. 6 a) 

(A. 6 b) 
(A. 6 c) 
(A. 6 d) 
(A. 6 e) 
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We are now in a position to solve the conservation equation for coupled matter in the 
Einstein frame. The Jordan frame stress tensor is conserved, as the matter it describes is 
uncoupled in this frame, so we can instantly write down: 


= 0 , 


(A.7) 


where is the covariant derivative metric compatible with g^ v . Given the transformations 
(A.6), (A.7) reduces to: 


p + a 


0 + 


3 

2 C 




(A.8) 


which, as long as w is constant, is exactly solvable: 

poc — (C' 1 / 2 a)- 3(1+u,/72) . (A.9) 

7 

That w be constant is not as restrictive a requirement as it sounds. In fact, as matter in 
the Jordan frame is uncoupled from the scalar, we expect w to be constant wherever it is 
in ACDM. For example, one can show that for any relativistic species (photons, massless 
neutrinos. . .), w = 1/3. 


B Frame transformations: perturbations 

Derivation of the transformations between perturbation variables of the two frames proceeds 
in the exact same way, though this time we will not be able to solve the equations exactly, as 
it can not be done for the uncoupled case. 

As discussed in section 3, <5 will denote matter density contrast: S := 5P the pressure 

perturbation, S(j) is the perturbation of the scalar held and 5 a general perturbation operator. 
Then, working in the Newtonian gauge to first order: 

ds 2 = a 2 (r) [-(1 + 2T)dr 2 + (1 - 2<S>)6 ij dx i dx j ] (B.l) 


which means that: 


ds 2 = Ccl 2 (t){ — (1 + 2A)j 2 dT 2 
+ 2(diB)'ydTdx l 
+ (1 — 2 E)5ijdx l dx^} 


where we recall the definitions of A, B, and E: 


. T SC Sj 

A = q,+ 2c + i 

1 \ S(j) 


B \7 7 M 

2 C ’ 


(B.2) 


(B.3) 

(B.4) 

(B.5) 


- 25 





where 


(B. 6 a) 

(B. 6 b) 


J7 

5C 


d r ZD' 

"7 c 

C'5(j) . 



X8(j) + SX 


Note that 4 1 — & = A —— As a consequence, if in the Einstein frame 4> — 4/ = 0, 
implying that the gravitational slip 7 = P/T = 1, the slip in the Jordan frame fj = E/A will, 
in general, depend on the coupling. 

Perturbations to the tensors given in (A.2b) and (A.lb) respectively gives: 

T» v = ( 8p + SP^v? + 2 (p + P)5u^u u) + SPg^ + P8g^ + PIP 1 ' (B.7a) 

V lv = (dp + 5P)u^u u + 2 (p + P)5u^u u) + 5Pg llu + PSg+ PIP 1 ', (B.7b) 


where we have denoted the fluid’s anisotropic stress IP"', which parameterizes higher moments 
of the fluid decomposition. And, just as before, by perturbing the map (A.3) we can compute 
the transformations between perturbed matter variables: 


8 

(p + P)~e 

8~P 

fp 1 ' 


s + - 
1 

c 2 

1 

C 

IP 1 ' 


2 5C J 7 

-- + 7 

(p + P)S+( 



7 ) 



(B. 8 ) 

(B.9) 

(B.10) 

(B.ll) 


These equations agree with [43]. 9 is the velocity divergence held, defined as 9 := di(ci5u l ), 
and quantifies discrepancies between the fluid’s velocity held and the underlying geodesic 
held. We have also defined for dark energy 9$ := Mathematically, every held permits a 
fluid description under a change of variables, and 9$ is the scalar’s velocity divergence. 
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